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In this paper the question about statistical properties of block-hierarchical random matrices is 
raised for the first time in connection with structural characteristics of random hierarchical networks 
obtained by mipmapping procedure. In particular, we compute numerically the spectral density of 
large random adjacency matrices defined by a hierarchy of the Bernoulli distributions {gi,q2,---} 
on matrix elements, where depends on hierarchy level 7 as = p~'^'^ {fi > 0). For the spectral 
density we clearly see the free-scale behavior. We show also that for the Gaussian distributions on 
matrix elements with zero mean and variances a.y = p""^^, the tail of the spectral density, pg{X), 
behaves as pg(A) ~ fgj- |^| _^ qq gj^^ < u < 1, while for u > 1 the power-law 

behavior is terminated. We also find that the vertex degree distribution of such hierarchical networks 
has a poly-scale fractal behavior extended to a very broad range of scales. 

PACS numbers: 05.40. a, 87.15.hg 



I. INTRODUCTION 

The determination of the structural organization of the complex system with statistical disorder refers often to 
a "network paradigm" , implying the investigation of topological and statistical characteristics of the contact map 
between the system elements. Analytical approaches to this task commonly use the well-defined relationship between 
the spectral properties of the adjacency (connectivity) matrix, which codes the contacts, and a particular topological 
structure of the network. In the case of uncorrelatcd random networks (graphs) known as Erdos-Rcnyi (ER) graphs 
there exists a bunch of methods for analytical study of spectra of the adjacency matrix - from mathematically 
rigorous to less rigorous. The last ones are often borrowed from the statistical mechanics of disordered systems. 
Among the most used are the replica or the Flory mean-field 0] approaches. 

In the last decade, an explosive growth of computational power has allowed to accumulate a large body of data 
on statistical and topological characteristics of real networks of diverse nature. The "network paradigm" is currently 
spread from protein folding, intermolecular contacts in biopolymers, genetic maps and cell metabolism, up to natural 
networks (including world wide web, various ecological, social, financial and economic entities, etc.), statistical data 
analysis, and even to Bose-Einstein condensation. Many of these topics are reviewed in In majority of cases 
it turned out that the statistical characteristics of real network topology, as well as the spectral properties of the 
adjacency matrices, essentially differ from the characteristics of random ER-graphs. For example, it has been found 
that the probability distributions of typical topological characteristics (e.g. the vertex degree distribution, or the 
clustering coefficient) have the power-law tails for many real networks, in contrast to the exponential tails for ER- 
graphs. Thus, the networks of such topology were stood out for a class of the "scale-free" networks just because 
of the power-law behavior of some observables. In practice, the network is often referred to belong to the "scale- 
free" class if the vertex degree distribution (i.e. the distribution of nearest-neighboring contacts over the network 
vertices) has the power-law tail. However in general setting dealing with the analysis of the spectral density of the 
adjacency matrix we shall understand under the "scale-free" -behavior the existence of the power-law tail in the 
spectral density. The topological characteristics of these networks, such as, for example, the vertex degree, have 
extremely wide distributions exceeding the corresponding distributions for ER-graphs by orders of magnitude. Such 
anomalously wide distributions we shall call "poly-scaled" to distinct them out of the "scale-free" behavior. 

The main idea of our work consists in the construction of networks with block-hierarchical adjacency matrices 
resulting from an appropriate randomization of the standard Parisi matrix - one of the key objects in the theory of 
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spin glasses (see @ for example). In other words, we study below the random block-hierarchical (RBH) networks. 
By the behavior of the spectral density of adjacency matrices, the RBH networks fall into the "scale-free" class, but 
in majority of cases, they have "poly-scaled" distribution of the vertex degree. 

It would not be an exaggeration to say that the block-hierarchical "ordering" is rather typical than exceptional for 
many complex systems which are both random and multi-scaled. The examples of such networks can be easily find 
in different areas of mathematics, physics and biology: from chaotic maps in Hamiltonian systems 0, 0] to condensed 
(globular) structures of polyelectrolyte chains [1| and hierarchical organization of biopolymers Q . We anticipate our 
consideration of spectral properties of the RBH adjacency matrices with a generic example of possible realization of a 
block-hierarchical contact map in a globular phase of unknotted ring polymer molecules with topological interactions. 

It is known that the non-phantomness of a polymer chain causes two types of interactions: i) volume interactions 
vanishing for infinitely thin chains, and ii) topological interactions, which present even for chains of zero thickness. For 
sufficiently high temperatures, a polymer molecule strongly fluctuates without reliable thermodynamic state called 
a coil state. However for temperatures below some critical value, 9, the poly meric chain exhibits a dense weakly 
fluctuating globular (drop-like) structure [lH • In classical works [13, devoted to the coil-to-globule phase 
transition without topological constraints, it has been shown that, for T < 6, the globular state can be described by 
using just only two- and three-body interaction constants: B = ^^-^-^ < and C = const > (see [H, [l^). The 
approach developed in [lol . [llj | is regarded as the basic one in modern statistical theory of collapsed polymeric state. 

The topological constraints in the globular phase of an unknotted macromoleculc act the part of repulsive interac- 
tions. For the temperatures below 0-point (i.e. in a poor solvent), against the temperature and the energy of volume 
interactions there exists certain scale of the chain length, g* , such that the chains longer than g* collapse. Taking 
an enough long chain, we can define these g*-unit parts as new "block monomers" (or "folds of minimal scale"). In 
the Figll^ they are denoted as the 1st level folds. Sufficiently long parts of the chain with several folds of minimal 
scale should again "collapse in itself", i.e. they should form the 2nd level folds, if other chain parts do not interfere 
with it. The chain of such new sub-blocks of the 2nd level folds collapses again forming the 3rd level folds, and so 
on... This block-hierarchical folding is completed when the initial chain units are united into one fold of the largest 
scale. Three first consecutive steps of such a process are shown in Fig[T^. Note that the line representing the chain 
folded by this way resembles the 3D-analoguc of the well known self-similar Peano curve. The specific feature of the 
crumpled globule consists in the fact that different chain parts are not entangled with each others, completely fill the 
allowed volume of space and are "collapsed in themselves" starting from the characteristic scale g* . 

The scale of the fold, shown by hues of gray in Fig[T^, can be considered as a cutoff for the interaction distance 
between (7*~block monomers in the current fold. The values t^"^ can be thought as interaction constants between 
these (;*~block monomers in the 7-level fold and, thus, they constitute the contact map T depicted in Fig[T}j. 

In [13] and later, more rigorously, in [31 it has been argued that the absence of knots in a densely packed polymer 
ring causes a very peculiar fractal structure of the chain trajectory, strongly affecting all thermodynamic properties 
of the macromolecule in the globular phase. The corresponding structure of a collapsed unknotted polymer ring was 
called a crumpled globule. The chain trajectory in the crumpled globule densely fills the volume such that all part 
of the chain become segregated from each other in a broad region of scales. This model has been used later (see, for 
instance, [IB, [l3|) to describe self-similar hierarchical organization in some biopolymers, like DNA and chromatin. 

Below we discuss some topological properties of block-hierarchical random networks. On the basis of obtained 
results we propose the new way of building of hierarchical networks with scale-free and poly-scaled properties in a 
memoryless locally uniform way. This method does not demand the control of the current state of the network. 



II. "HEAVY TAILS" IN THE SPECTRAL DENSITY OF RANDOM BLOCK-HIERARCHICAL 

ADJACENCY MATRICES 

Let us start with a description of generic procedure of the RBH-network construction. Taking N points as potential 
vertices of our forthcoming network, we raise a hierarchical network by connecting the vertices by edges in a specific 
way. In the outset, we introduce an ensemble of the N x N adjacency matrices, T; each of them encoding the 
edges between connected vertices in a network realization. Namely, the element Tij of T is 1 if the vertices i and 
j are directly connected, otherwise Tij = 0. Wc consider the adjacency matrix in very peculiar form of a p-adic 
translation-noninvariant Parisi matrix. This matrix is shown in Fig[T]3 for p = 2. Obviously, Tij = Tj^i and Ti^i = 0. 
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(a) (b) 

Figure 1: (a) Three subsequent stages of the construction of hierarchical contact map for the "crumpled globule" (see the text 
for details); (b) Block-hierarchical p-adic Parisi matrix T (p = 2). 



All matrix elements, tx^\ are the Bernoulli distributed random variables: 

I 1 with the probability 
I with the probability I — 

where 7 counts the hierarchy levels (1 < 7 < 7max = T) and n enumerates different blocks corresponding to a given 
hierarchy level 7 (see FigUlD). Note that the probability q^ does not depend on n. The full ensemble oi N x N 
matrices T, where N = , is completely determined by the set of probabilities, {Q} = {gi, 52, 9r}- Thus, the 
elements Tij, being the random variables, are hierarchically organized in probabilities. Below we consider the set of 
probabilities, {Q}, with q^ = p^'''' (^ > 0). 

The systematic study of statistical properties of ensembles of random graphs (networks) deals with the investigation 
of the spectral properties of a graph adjacency matrix [l^, [l^. Let Xi {1 < i < N) be the eigenvalue of the adjacency 
matrix. The spectral density of the ensemble of random symmetric adjacency matrices is defined in the standard way, 

1 ^ 

P(A) = ^EWA- A.)) (2) 
1=1 

where {■■■)^g^ g y denotes the averaging over the distributions of the matrix elements, t^\ 

Computing numerically the spectral density, p(A), of networks with block-hierarchical adjacency matrices, we found 
that the tails of the spectral density p(A) follow a power-law asymptotic behavior p(A) ~ |A|^-*^ with the exponent 
X = x(m)- The sample plots of the spectral density p(A) are shown for N = 256, 2048 and yu = 0.2 in Figl2^ in semi-log 
coordinates. The corresponding log-log plot of the left- and right-hand tails of the spectral density for N = 256 is 
drawn in Figl^ja. It is interesting to note that the right-hand tail of p(A), while demonstrating the same behavior, is 
worse averaging. 

Some analytic arguments supporting the found power-law behavior of the spectral statistics of the RBH- graphs 
can be brought by means of the spectral density, Pg(A), of the Gaussian ensemble of the Parisi matrices T, where the 
Bernoulli distributions on matrix elements of T are replaced by the Gaussian distributions with zero means and a set 
of variances {S} = {tTi,cr2, ...,o'r}, cr-^ = p^'^'^ {v > 0). To get the spectral density pg(A), note that the eigenvalues 
of the standard (i.e. translation-invariant) Parisi matrix with Tt^i = can be expressed in terms of matrix elements 
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(a) (b) 

Figure 2: Spectral density p(A) for /i — 0.2: (a) semi-log plot of the full distribution for N = 256 (solid line) and A'^ = 2048 
(dashed line), (b) the left- and right-hand tails of p(A) for A*" — 256 in log-log coordinates. 



4"^ = as follows (0): 

7 

A^=p^t^-(l-p-i)^p'''iy (7 = l,...,r) (3) 

7' = 1 

The eigenvalue is p^~^ times degenerated (7 = l,...,r). In addition there is one extra eigenvalue Xq = 
~Y^^=iP^~''^f For the translation-noninvariant Parisi matrix we generalize ([3]) in the way similar to the one 
used for block-hierarchical kinetic matrices (see [2l|). Remind that for eigenvalues of kinetic matrix, one has 

r 

A7,„ = -p^t(")-(l-p-i) J2 p'''*y'^ (4) 

7'=7+l 

S 

/ /\ 

i.e. the eigenvalue A^.„ can be expressed via a linear combination (E in Q) of weighted matrix elements ty coming 
from the vertices (7',^') along a unique path on the p-adic Cayley tree from the vertex (7' =7,71' = n) towards 
the root vertex (7' = F, n' = 1) - see [22| for more details. In the framework of same geometrical interpretation, the 
eigenvalue A^.„ of the block-hierarchical matrix with T^ ^ = (shown in Fig[T|3) seems to be the linear combination 

of weighted matrix elements along a path on the p-adic Cayley tree from the bottom level 7' = 1 to the 
vertex (7' = 7,n' = n). However this construction has an ambiguity since the corresponding path is not imiquely 
defined by the pair (7,71). So, rigorously speaking, the eigenvalues of particular realization of the random block- 
hierarchical adjacency matrix cannot be parameterized by the pairs (7,n), and we cannot write an exact expression 
for the eigenvalues in the form of expression ([3]). However for computation of spectral density, pq{X), we use below a 
posteriori self-averaging arguments, which make our consideration self-consistent. Moreover, the extensive numerical 
simulations confirm our analytic prediction of pcW for |A| 3> 1 in the interval < < 1. Thus, formally extending 
^ to the case of A-y_„ of block-hierarchical matrix T, we replace in ([3]) the first term by t^"^ and the sum - by 
X^y^i P'^ ^7" ^ where the summation runs now along the paths on the Cayley tree from the hierarchical level 7' = 1 to 
the vertex (7,?^) located on the hierarchical level 7' = 7 (compare to (|31)). Supposing the distribution of the matrix 
elements t^'' to be Gaussian, 



P(t(")) = ^exp^-i^) (5) 



M\2 ' 



and using for A^.„ the guessed expression as a linear combination of matrix elements, we end up with the following 
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equation for the spectral density, pq{X): 



1 ^ _ 1 



7,n 



7=1 \ 1^' 



cxp 



where for = p '^'^ we have: 



7' = 1 



P — P'" 



For p = 2 we can rewrite pcW in ® for F — > oo as follows 



oo 

Pg(A)c^^5]2-(2-)7cxp 



7=1 



4-4"^ 



22(1-1^)7 _ 1 



Taking into account that 



7=1 



(6) 



(7) 



(8) 



(9) 



and substituting ci = 2 ~ v and C2 = 1 — i^, we arrive at the following asymptotic form for the spectral density at 
|A| » 1: 



Pg{X) ^ IM-^^"'^ (0<i.<l) 



where 



2-v 
l-v 



(10) 
(11) 



The arguments supporting our derivation of the expression ([7]) for the spectral density pg(A) are as follows. First of 

all, note that ©-([7]) become exact if we skip the dependence on n in the matrix elements t^"^ and, hence, restore the 
translational invariance in the block-hierarchical matrix T. Secondly, we found in the extensive numeric simulations 
summarized in Fig[3]that indeed the conjectured behavior (|10p - pip actually holds for translation noninvariant Parisi 
matrices. In Fig[3]we have plotted the tails of the spectral density pcW for Gaussian ensemble of block-hierarchical 
matrices for TV = 256. The solid and dot-dashed lines have the slopes — C(i'), where ^(i/ = 0) = 2 (solid line) and 
^(i/ = 0.8) = 6 (dashed line). The scatter graphs from top to bottom correspond to ly = 0.0,0.2,0.4,0.6,0.8, 1.0. 

The fact that pU]) gives right asymptotic behavior occurs apparently due to an effective self-averaging of the sum 
of matrix elements along each particular path on a Cayley tree for the distribution {S} = {ai, (72, ■■■}■ One sees 
that in the sum in ([7]) for a-y = p~'^'^ the lower limit of the summation can be shifted from 7' = 1 to 7' ^ —00. 
Asymptotically the result for it^ will remain unchanged if < < 1. Such an extension of summation means that 
the computation of the spectral density (and, in particular, of A^.„) involves the summation along the infinite paths 
running from ^00 to the hierarchical level 7. We expect that for = p~'^^ (0 < < 1) due to the convergence of the 

sum J2y=~ooP'^'^ '^Y' eigenvalue A^_„ does not depend on each particular path on a Cayley tree and hence, the 
eigenvalue does not depend on the index n. Once this point of view is accepted, we return to an effective translation 
invariant block-hierarchical matrix for which ([5]) is exact. 

One can see from ([TO)) -pT |) and Figl3]that for > 1 the power-law behavior of the spectral density pg(A) terminates. 
This termination deserves special attention. Indeed, for v > 1 wc cannot extend the lower limit of summation over 
7' to —00 since the corresponding sum diverges. Hence, the contribution to the eigenvalues (and, therefore, to the 
spectral density) strongly depend on the particular configuration of the path. In this case the self-averaging, suggested 
in our computation, is invalid anymore and we cannot say anything about the behavior of the spectral density in the 
region v >1. The opposite case v <0 deserves special attention because formally the sum ([7]) converges for any v <{) 
and our arguments about self-averaging seem to work. However the details of the analysis of this case is beyond the 
scope of the current work and will be discussed elsewhere. 

Our approach to the construction of the networks with the scale-free behavior in the spectral density is not unique. 
As it is shown in the recent work [23j . the procedure of dividing Gaussian matrices by a random variable, as well as 
the same procedure applied to random graphs, leads to the spectral density interpolating between the Erdos-Rcnyi 
and the scale-free models. 
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Figure 3: Spectral density for symmetric Gaussian translationary noninvariant block-hierarchical matrices. The solid and 
dot-dashed lines have the slopes ^{v — 0) — 2 and = 0.8) = 6; scatter graphs; v = Q (+), y = 0.2 (0), v = 0.4 (t), v = 0.6 
(A), j.^0.8 {•),u^l.O (■). 



III. HIERARCHICAL "GROWING" OF POLY-SCALED NETWORKS 



One can roughly distinguish two methods of construction the scale-free networks. The first method has been 
developed mainly for illustrative purposes and deals with the hierarchical construction [13] of deterministic scale-free 
graphs with predetermined fractal properties such as, say, vertex degree distribution. Since in this case the obtained 
graphs are deterministic, it is senseless to talk about any statistics of their spectra. The second method deals with 
the variants of iterative "preferential attachment" construction (H, where new nodes are added to a vertex of the 
network with the probability depending on already existing vertex degree (the number of one-step connections to 
other vertices). Almost all known statistical characteristics of scale-free networks, including the spectral density of 
adjacency matrix, are obtained for the networks constructed using this method. Typically, the spectral density of 
the ensemble of scale-free networks designed by the preferential attachment method has a triangle-like shape in the 
"bulk" part with power-law tails. 

Note that the construction of scale-free networks by preferential attachment method is based on locally nonuniform 
incremental growth with unlimited evolutionary memory. In contrast to this, we propose below another physically 
motivated approach to the construction of the random networks with power-law spectral density by a "parallel" (i.e. 
non-stepwise) and uniform procedure. 

The peculiarity of our construction consists in the following. We build clusters of edges with hierarchically organized 
probabilities, while the typical procedure consists in hierarchical grouping of vertices. Our procedure does not impose 
any additional metric structure on the graph and leaves the graph (network) purely topological. To the contrary, 
grouping of vertices imposes a metric structure on the graph (network) since such a grouping operates usually with 
the notion of "close" (or "distant") vertices. 

To shed light on spectral properties of random hierarchical graphs we exploit the link between the Ran- 
dom Matrix Theory (RMT) and the Graph Theory (GT). It is known that the spectral density p{X) = 
{i3i,g2,- -,<?r} ensemble of adjacency matrices is directly related to the topological structure of 

the corresponding network since the value 

1 /■ 1 ^ 

i—l 

defines (up to the factor N) the average number of /c-stcp loops in the network (sec, for example, [l8|). 

In particular, in [26j it has been shown that in the thermodynamic limit the spectral properties of random Erdos- 
Renyi graphs [27| coincide with the spectral properties of random real symmetric matrices. This result is one of 
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the benchmarks in our consideration. The elements of adjacency matrix, Aij, of Erdos-Renyi graph are BernouUi- 
distributed random variables: A^j = 1 or with probabilities q and I ~ q correspondingly. For the ensemble of 
ER-random graphs the density, paW, of eigenvalues of adjacency matrices A can be analytically computed in the 
thermodynamic limit iV — > oo and has for some q a celebrated Wigner-Dyson semicircle law known for ensembles of 
Gaussian matrices [1^. Namely, the following statement is proved [26|. Let _B be a real symmetric N x N matrix 
with independently distributed entries Bij from, say, Gaussian distribution P{Bi,j) with {Bi,j) ~ and (Bij) = a". 
Then the spectral density. pb{X), of the ensemble of matrices B converges in the limit cx) to the semicircle 

distribution 



^y/ANcr^ - \^ if |A| < ViNa^ 



Pb{\) = <! ■ ■ (13) 

[o if|A|>V4]V^ 

If cr^ = q{l — q), then p/i(A) = Pb{^) for N oo, i.e. the spectral densities of ensembles of random ER-graphs and 
of Gaussian symmetric matrices coincide in the thermodynamic limit. Nevertheless, such a coincidence of spectral 
densities for random graphs and random Gaussian matrices should not be understood in a literal sense: some spectral 
properties of random ER-graphs and random matrices are different p^ . [29| . For example, since for the adjacency 
matrix A of random ER-graphs one has (Aij) = q, then the corresponding largest eigenvalue, Ai, grows linearly with 
the system size, N, i.e. Ai = Nq, meaning that the semicircular distribution for random graphs is valid only for the 
matrix A — (A) . Also, the tails of spectral distributions near the spectrum edges are different for random graphs 
and random matrices. Nevertheless, as a first approximation, the random Gaussian matrices could serve as a very 
natural benchmark for the corresponding statistical analysis. The comparison of spectral properties of hierarchical 
and random Erdos-Renyi graphs is demonstrated in Figdl The FigH^ shows the semi-log plot of the spectral densities 
for: (i) random hierarchical graphs with LCP adjacency matrices for N = 256 and p = 0.2 (solid line), and (ii) random 
Erdos-Renyi graphs for N = 256 and p = 0.2 (dotted hue), for N ~ 256 and -p = 0.02 (dashed line). In Figl4)D we 
have redrawn the central part of Figd^ in the linear scale. 




Figure 4: Comparison of spectral densities of random block-hierarchical matrix and Erdos-Renyi (ER) random graphs: (a) in 
semi-logarithmic coordinates in the region —0.5 < A/A^ < 0.5; and (b) in linear coordinates in the region —0.05 < A/A^ < 0.05. 



The numerical results for the probability distributions P{Mk) {k = 2, 3) of the number of fc-step loops in random 
hierarchical graphs for N = 256 and p — 0.2 are compared in FigOwith the distributions of loops on Erdos-Renyi 
random graphs for N = 256 and p = 0.2. Recall that P{M2) defines the probability to have in the finite graph the 
average connection degree equal to M2. 

One sees that the distribution functions P{M2) and P{M3) for our hierarchical random graphs arc much broader 
than the corresponding distributions for Erdos-Renyi graphs with the same number of vertices. Hence, the topological 
structure of random hierarchical graphs is much more "flexible" than that of random ER-graphs. This is consistent 
with the found behavior for the spectral density: the distribution function p(A) has "heavy" tails and decays much 
slower than the that for random ER-graphs. According to the behavior of the distribution functions P{M2) and 
P{M-i) it is naturally to call our random block-hierarchical graph the "poly-scaled". 
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Figure 5: Distribution of the number of 2- and 3-step loops in random hierarchical and Erdos-Renyi graphs. 



The method of generating functions allows us to compute easily the vertex degree distribution in the ensemble of 
random block-hierarchical graphs directly by their adjacency matrices. To do that let us consider any (for example, 
the first) row in the adjacency matrix T (see FiglT})). The total number of links of the first graph vertex to other 

vertices is defined by the number of matrix elements, t^^^ (7 = l,...,r), having in the first row nonzero value (i.e. 
taking the value "1"). Thus, the distribution of the number of connections (i.e. vertex degree distribution), V{m), is 
the probability of the fact that the sum of matrix elements in the first row is exactly equal to m under the condition 
that the matrix elements are grouped in the hierarchical blocks and have the binomial distributions {91, 92, ■■■,<lr} 
as it is defined in ([T]). Finally we arrive at the following expression for the degree distribution V{m) (in the sake of 

simplicity we have denoted q-^it^^) = g^(t-y)): 



V{m)= J2 



{ti...tr} 




where the binomial distributions q-y{t-y) have the form 

q^m)^p-^''S,^^t,)^, + (1 -p-'^7)5,^(,^),o 

and A(...) is the Kronecker symbol: 



(14) 



(15) 



A(a:) ^ — (tdz 



1 if x = 
if x 7^ 



Substituting (|15p and in after elementary transformations we get 



where 



■^'^^ 7=1 



(16) 



(17) 



(18) 



For not too large values of F the distribution V{m) can be analyzed numerically. Using the fact that the function 
W(z,7) is a polynomial of z, let us represent Vijn) as: 



Vim) 



1 



jm+l 



(m+ 1)! 



(19) 



z=0 



In Figin^ we have depicted the family of curves P{m) for F = 16 and ^ = 0.1; 1.0. For comparison, in FiglSJ; we have 
plotted the distribution 7^(m) for /i = 1.0, F = 16 as well as the binomial distribution 7'er("^) = (7™(1 — g)^"™ 



On scale-free and poly-scale behaviors of random hierarchical network 9 




Figure 6: a) The family of distributions Vijn) for F = 16 at = 0.1; 1.0; b) Distribution Vim) in the double-logarithmic scale 
for ^ = 0.1 and F = 16; c) Comparison of 'P{m) for the random block-hierarchical graph — 1.0, F = 16) and PerC'ti) for 
the Erdos-Renyi graph {q = 0.1; 0.5, A*' = 2^**). 



for the standard Erdos-Renyi graph with the number of vertices N = 2^ = 2^^ and for two values q = 0.1; 0.5. 
One can see on FiglSJ; how wider is the distribution V(rn) with respect of the corresponding distribution ■Per(to) for 
ER-graphs. 

The fractal structure of the distribution T'{m) for random hierarchical networks demonstrated in Figl6^,b, ap- 
parently is deeply linked to the invariant multifractal measures appearing in chaotic Hamiltonian systems in 
connection with the problems of "number-theoretical chaos" - see, for example, [sO]. Actually, the condition 

A ^X]7=i P'^i-y ~ ™ equation (jl4p for F — > oo is nothing else as the binary expansion of the number to: 
TO = ti 2° -|-t2 2^ -l-ts 2^ + ... 2"' + where the coefficients take values 1 or with corresponding probabilities 

q.y{t^) defined in (jlSp . Let us note that the similar expansion of the form X]fc°=i ^kU~'', where u > 1 and £k = ±1 
(with equal probabilities i independent of k) is known in the literature as the "singular Erdos measure" 
The observed fractal structure of the distribution function V{m) appears due to the effects of the incommensurability 
of number-theoretic origin: some binary expansions (of the number to) with random coefficients have relatively high 
probabilities to appear in the adjacency matrix, while other binary expansions have much less possibilities. 



IV. CONCLUSION 



First of all, let us note that our consideration of spectral properties of random hierarchical graphs is far from being 
complete and many other properties are of interest (for example, eigenvectors, inverse participation ratio, etc. — see, 
for instance [Toj). However even in this preliminary investigation we would like to emphasize the crucial difference 
between random hierarchical and random Erdos-Renyi graphs. 

Roughly speaking, there are two generic ways of the scale-free network construction. The similarities and differences 
of these two approaches we would like to emphasize below again. 

• The first way, widely discussed in the literature, deals basically with the "preferential attachment" procedure, 
where the network is raised by the essentially non-Markovian (in the increments) evolution process with unlimited 
memory. The evolutionary process of such kind can be viewed, to some extent, as a "generalized Brownian 
motion" . The realization of such a procedure demands the monitoring of the whole network structure on each 
step, because the appearance of new links depends on the current degree of graph vertices. From this point 
of view the sequential construction of the corresponding network can be tentatively denoted as a "nonlinear 
evolution" . 

• The second way, discussed in the present work, exploits essentially different mechanism of the scale-free network 
formation. The hierarchical organization of probabilities of links in the topological network constitutes the basic 
idea of our construction. Specifically, we construct the networks with the scale-free eigenvalues distribution 
of the adjacency matrices, where the last originates from an appropriate randomization of the standard Parisi 
matrix - one of the key objects in the theory of spin glasses. Our method allows one to build a hierarchical 
network in a memoryless locally uniform way. The application of the elements of the p-adic analysis [s^ permit 
us to analyze the basic spectral properties of ensembles of randomized Parisi-type adjacency matrices. 
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To end up, let us emphasize that there are two important features of hierarchical networks constructed in our paper. 
First of all, any sub-graph belonging to particular hierarchy level, 7, is just a random Erdos-Renyi graph because 
the formation of clusters of bonds on each hierarchy level is entirely uncorrelated. Secondly, the random sub-graphs, 
associated with different hierarchy levels of the network, can be different, so the network in a whole can be essentially 
nonuniform. Nevertheless as we have seen, the "mipmapping" construction with different sets of parameters (i.e. the 
hierarchical embeddings of sub-graphs corresponding to the different hierarchical levels) leads under some conditions 
to the scale-free behavior in the spectral density of adjacency matrix and to the poly-scale (and even to the fractal) 
behavior in various topological characteristics of the graph. This observation is rather unexpected since in our case 
the scale-free behavior is reached by essentially Markovian and mcmorylcss procedure. 

In physics the random graphs of such a hierarchical genesis can be encountered among the scale-free networks 
whose natural origin are associated with low-correlated random events carried out under short evolutionary memory. 
In particular, the networks of hierarchical genesis may by interesting for the prebiology or the earliest biology. 
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